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Abstract. We have mapped the large-scale structure of the Serpens cloud core using moderately optically thick 
(^^CO(l-O) and CS(2-1)) and optically thin tracers (C^**O(l-0), C^^S(2-1), and N2H+(l-0)), using the 16-element 
focal plane array operating at a wavelength of 3mm at the Five College Radio Astronomy Observatory. Our main 
goal was to study the large-scale distribution of the molecular gas in the Serpens region and to understand its 
relation with the denser gas in the cloud cores, previously studied at high angular resolution. All our molecular 
tracers show two main gas condensations, or sub-clumps, roughly corresponding to the North- West and South- 
East clusters of submillimeter continuum sources. We also carried out a kinematical study of the Serpens cloud. 
The i^CO and 0^*0(1-0) maps of the centroid velocity show an increasing, smooth gradient in velocity from East 
to West, which we think may be caused by a global rotation of the Serpens molecular cloud whose rotation axis 
is roughly aligned in the SN direction. Although it appears that the cloud angular momentum is not sufficient for 
being dynamically important in the global evolution of the cluster, the fact that the observed molecular outflows 
are roughly aligned with it may suggest a link between the large-scale angular momentum and the circumstellar 
disks around individual protostars in the cluster. We also used the normalized centroid velocity difference as an 
infall indicator. We find two large regions of the map, approximately coincident with the SE and NW sub-clumps, 
which are undergoing an infalling motion. Although our evidence is not conclusive, our data appear to be in 
qualitative agreement with the expectation of a slow contraction followed by a rapid and highly efficient star 
formation phase in localized high density regions. 
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1. Introduction 

It is now well established that a substantial fraction of 
the young pre-main sequence stars in molecular clouds 
are found in compact groups (clusters) rather than being 
uniformly distributed throughout the clouds (e.g. Gomez 
et al. |I993t Testi et al. |1999| Carpenter |00§). It is thus 
clear that the formation and early evolution of stellar clus- 
ters hold the key for the understanding of the galactic 
disk stellar population, such as the distribution of stellar 
masses at birth, or initial mass function (IMF), as well 
as the kinematical and multiplicity properties (Clarke et 

pooo| 



al 



Kroupa 



1995 



Kroupa et al. |2001| ). While there 
is a general consensus that the latter two properties are 
linked to the clusters dynamical evolution, the nature of 
the IMF is still debated. Observations show that the emer- 
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gent mass distribution of young stellar clusters is con- 



sistent with the field st ars IM F (Palla & Stabl er |1999 
Hillenbrand k Carpenter [2000[ Meyer et al. ^OOOD , imply- 



ing that the stellar mass distribution must be linked with 
the cluster formation process itself. Three main classes of 
models have been considered: those which link the final 
stellar masses to (i) the structure or fragmentation of the 



parent molecular cloud (e.g. Elmegreen 1997; Myers 2000 
Padoan & Nordlund |2002| ) ; or (ii) the feedback of the pro- 



tostellar accretion process (Adams & Fatuzzo 1996); or, 
finally, (iii) competitive accretion between already formed 
(proto-)stars in clusters (Bonnell et al. 1997; 2001). 



The progressive improvement in millimeter wavelength 
techniques, in particular the development of bolometer 
arrays for large single dish telescopes and the upgrade 
of millimeter wave arrays, have provided a new wealth 
of observational data to constrain theories (Andre & 
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Motte |2000| Testi & Sargent pOOOj ) . Millimeter continuum 
surveys of the cluster forming clouds in Serpens (Testi & 
Sargent |l99|, hereafter TS98), p-Oph (Motte et al 



1998 



Johstone et al. gOOO|), NGC 133 3 (Sa ndeU & Knee EOOlh, 
and Orion B (Johnstone et al. ^00 1| Motte et al. ^OOlD 



have shown that the mass distribution of the prestellar 
cores, progenitors of individual stellar systems, is consis- 
tent with the field stars IMF and thus significantly steeper 
that the mass function of the larger gaseous clumps in 
molecular clouds (e.g. Williams et al. pOOC ). These obser- 
vations suggest that the self-similarity of molecular clouds 
breaks down at the cores level, and the fragmentation or 
the small-scale structure of the clouds is the likely mech- 
anism responsible for the stellar IMF. 

In a first attempt to link the small-scale structure 
of the prestellar cores and young stellar objects in clus- 
ters with the kinematics and large-scale distribution of 
the molecular gas, Testi et al. ( |2000| ; hereafter TSOO) 
combined the high resolution 3 mm continuum survey of 
TS98 with new interferometric and single dish molecu- 
lar line surveys of the Serpens core. The high resolution 
CS(2-1) observations showed high velocity molecular gas 
associated with outfiows emanating from several of the 
protocluster members, while the single dish N2lI+(l-0) 
revealed the presence of two main gaseous clumps, well 
separated in velocity, in which the TS98 sources are em- 
bedded. The combination of the new high and low reso- 
lution millimeter- wave data with previously published ob- 
servations showed that the protocluster is structured in 
sub-clusters, each internally coherent, well separated both 
spatially and kinematically. The star formation efficiency 
and proto-stellar density within each core is a factor of 
~10 higher than the average values. The fact that the 
molecular gas within each sub-cluster is kinematically co- 
herent suggests that clusters are assembled through frag- 
mentation of the molecular clump in coherent sub-clumps 
which produce independent sub-clusters (TSOO). It is not 
clear, however, whether the sub-clusters in Serpens will 
eventually merge into a final cluster, such as other young 
clusters (e.g. the Orion Nebula Cluster) which show little 
evidence for sub-structures (Clarke et al. 200C| ). 

Similar results on the spatial segregation of pre-stellar 
cores have also been found in the p-Oph star forming re- 
gion (Johstone et al. 2000), but unfortunately the essen- 
tial kinematical information is missing in this region. In 
both regions, the prestellar cores and protostars densities 
are too low for dynamical interactions to be efficient over 
the timescale in which the bulk of the mass of the final 
stars is assembled. Dynamical interaction and competitive 
accretion can possibly play a role for the final tuning of 
stellar masses but cannot play a major role in determining 
the global shape of the stellar IMF in these regions. Note 
however that both these regions are at least an order of 
magnitude less dense than the typical clusters associated 
with massive stars (such as the Orion Nebula Cluster). 

Given for granted that, at least in Serpens and p-Oph, 
the stellar IMF appears to be linked to the individual 
small-scale cloud fragments, each progenitor of a single 



stellar system, the question remains on how the fragmen- 
tation and collapse processes occur. There are currently 
two opposite views of these processes: a "slow" and a 
"fast" mode. In the first mode molecular clouds may con- 
tract as a whole in a slow fashion evolving through a series 
of quasi-equilibrium conditions; star formation is triggered 
in localized regions of the cloud when the density exceeds 
the threshold for gravitational collapse. Star formation 
in the cloud as a whole may be regulated by external 
agents, such as in the photoionization-regulated view of 
McKee (1989), or it may progressively accelerate until the 
feedback of newly formed stars disperse the molecu lar ma- 
terial, halting star formation (Palla & Stabler 2000). In the 
"fast" mode, since turbulence is expected to decay o n ver y 
short timescales within molecular clouds (Mac Low 1999 ), 
there is no force to balance the clouds against gravita- 
tional collapse and star formation occurs globally on a few 
free-fall times (Elmegreen 2000 ; Hartmann et al. 2001 ). 
From the observational point of view, the problem can be 
approached by studying the large-scale kinematical and 
physical conditions in star forming molecular clouds and 
compare them with position, mass and age distributions 
of forming stars (Palla & Stabler 2000). 

The Serpens cloud core appears to be an ideal place 
to test the predictions of the various models. At a dis- 
tance of ^310 pc (de Lara et al. 1991), the Serpens cloud 
core and protocluster have been the subject of numer- 
ous detailed observational studies in recent years, at in- 
frared and (sub-) millimeter wavelengths. Among the most 
recent are: Davis et al. 1999: Giovannetti et al. 1998; 
Hodapp |1999[ Hog erheijde et al. |1999|; Kaas |1999| ; Kaas 
& Bontemps 200lt M cMuUin et al. [1994| |200C| ; Williams 
& Myers |1999| ; ^OOOt in addition to TS98 and TSOO. The 
reasons for these numerous observational efforts lie in its 
compactness, richness and relative proximity to the Sun, 
which make it possible to obtain and combine high resolu- 
tion observations of the relatively large number of (proto- 
)stellar sources within the (proto-)cluster and large-scale 
observations in various tracers of the entire cloud core. 
Regions further away from the Sun cannot be adequately 
resolved with single dish telescopes and cannot be ob- 
served at the required sensitivity with current millime- 
ter arrays, whereas nearby regions require a substantially 
larger observing time to be fully mapped. 

Most of the previous single dish millimeter studies of 
the large-scale structure of the molecular clump, within 
which the (proto-)cluster is embedded, focused either on 
the study of molecular outflows by means of optically thick 
C0(2-l) observations (Davis et al. 1999), or on the chem- 
istry (McMulhn et al. [1994| |2000D , or present ed pa rtial 
maps of the cloud core (White, Casali & Eiroa 1995). In 
this paper we complement these and our previous high 
resolution studies (TS98 and TSOO) with observations 
of the molecular clump large-scale structure using com- 
plete maps in moderately optically thick (^'^CO(l-O) and 
CS(2-1)) and optically thin tracers (01^0(1-0), C^'^S{2- 
1), and N2H+(l-0)), taking advantage of the 16-element 
focal plane array operating at a wavelength of 3mm at the 
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Table 1. Molecular transitions observed and their rest 
frequencies. The third column lists the size of the map cen- 
tered about the position a = 18^27™20" and 5 = 1°12'30" 
(B1950.0) 



Line 


Frequency [GHz] 


Map size [arcmin] 


CS(2-1) 


97.980968 


13 X 10 


C3*S(2-1) 


96.412982 


5x5 


N2H+(l23-0l2) 


93.1734796 


13 X 10 


"CO(l-O) 


110.201353 


16 X 16 


Ci«O(l"0) 


109.782182 


12 X 14 



Five College Radio Astronomy Observatory^ (FCRAO). 
These new observations are used to derive refined values 
of the physical parameters of the cloud, and are combined 
with previous observations to study the kinematical and 
fragmentation properties of the cloud core. 

The paper is organized as follows: in Sect. ^ and || 
we describe the new observations, the morphology of the 
Serpens cloud in the observed tracers and derive the phys- 
ical conditions of the molecular gas; in Sect. I we discuss 
the cloud kinematical and fragmentation properties; fi- 
nally, in Sect. |5| we summarize the main results of this 
study, and we discuss the implications for cluster forma- 
tion and evolution theories. 

2. Observations 

The single dish observations of the Serpens cloud core 
were carried out in March, October and November 1999 
with the 13.7-m telescope of the FCRAO located in New 
Salem (U.S.A.), using the new SEQUOIA 32 element fo- 
cal plane array, although for the present observations only 
16 pixels were actually used. The observed lines and their 
rest frequencies are listed in Table |l|. The maps were cen- 
tered about the position a = 18^27^20'^ and 5 = 1°12'30" 
(B1950.0), and the FWHM of the FCRAO telescope var- 
ied from about 46" to about 53". The front-end receivers 
employed low noise InP MMIC-based amplifiers result- 
ing in a mean receiver temperature of 70 K (SSB) and 
system temperatures typically in the range 160 — 250 K. 
Our spectrometer was an autocorrelator with 24 KHz 
(0.077 km s^^) spectral resolution and 20 MHz band- 
width. 

The integration times used for the single-scan were 
typically 5 to 25 minutes in frequency switching mode, 
using a frequency throw of 8 MHz. The maps were carried 
out using I-beam sampling and the sizes of the covered 
regions in the different lines are also listed in Table |l|. 
The main beam efficiency used to convert the antenna 
temperature, T^, to main beam brightness temperature 
is ?7mb = 0.55. The final achieved sensitivities in the 
scale varied from about 0.05 K for C'^'^S to about 0.19 K 
for 13 CO. 




1 The Five College Radio Astronomy Observatory is oper- 
ated with support from the National Science Foundation and 
with permission of the Metropolitan District Commission 



Fig. 1. Maps of the ^^CO, C^^O and N2H+ integrated 
intensity in the velocity interval 4.0 to 12.0 km s^^. In the 
13 CO map the value of the first contour is 1.5 K km s~iand 
the contour interval is 1.0 K km s~^. In the C^^O and 
N2H+ maps the value of the first contour and the contour 
interval are both 0.4 K km s~^. The cross marks the refer- 
ence position (see Table and the filled triangles repre- 
sent the positions of the submillimeter continuum sources. 
The letters mark the positions of the main N2H+ cores as 
defined in TSOO. The "spur" is also visible, an extended, 
lower surface brightness narrow region, stretching from 
the NW to SE clumps and continuing beyond these (see 
also Fig. I). 



3. Results 

3.1. Morphology of the Serpens region 
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3.1.1. "CO(l-O) and Ci*^O(l-0) 



°C and BCuaA 850 f^n-. 

T- 



We performed wide-field imaging of tlie ''CO and 
C^^O(l-O) emission towards the Serpens region to under- 
stand the relation of the global distribution of the molec- 
ular gas with the denser gas in the cloud cores. The opti- 
cally thicker CO (1-0) line emission tends to pick up the 
more diffuse ambient gas, whereas the thinner transition 
of C^*O(l-0) shows the distribution of the gas inside the 
Serpens cloud. 



The 13CO(1-0), Ci*O(l-0) and N2H+(l-0) emission 
integrated over the velocity interval 4.0 to 12.0 km s""'^, 
shown in Fig. ^ trace the distribution of the column den- 
sity in the Serpens cloud (see also Sect. 3. 4.1] ). The ^■^CO 
emission is extended and encompasses both the NW and 
SE clusters of submillimeter continuum sources, with a 
peak of emission near the position of SMM4. We clearly 
detect 1^00(1-0) at the position of S MMl , where the 
BIMA observations of McMuUin et al. ( |2000| ) did not de- 
tect any emission. Compared to the C^*^O(l-0) map of 
McMuUin et al. (2000) our C^^O integrated intensity map 



better shows the clear demarcation between the NW and 
SE sub-clusters, which are also separated in velocity (see 



Sect. 3.3 and TSOO). 



The NW clump is clearly surrounding all continuum 
sources: S68N (McMuUin et al. |1994| ) and SMMl in par- 
ticular (also seen at the position of the N2H+ core C) are 
very close to two C^^O emission peaks. The SE cluster 
of submillimeter continuum sources lies essentially in the 
region between the two C^^O peaks, and is aligned with 
the C^^O lower surface brightness region extending to the 
NE. On the other hand, sources SMM2 and SMMll are 
embedded in the SE gas clump, and the emission is elon- 
gated in the NS direction along the "spur" , an extended, 
lower surface brightness narrow region, stretching from 
the NW to SE clumps and continuing beyond the se, ob - 
served by TSOO. As observed by McMuUin et al. (|2000|) , 
in the case of the C^^O emission there is a lack of corre- 
spondence with the submillimeter peaks, especially in the 
SE sub-clump. However, the C^*0 and submillimeter dust 
continuum emission both trace the large-scale distribution 
of column density, as shown in Fig. ^. A similar structure 
to C^^O can be seen in N2H+ (TSOO), although the peaks 
of emission show a different position in the SE clump. 

Fig. H shows that the Spur emission has a striking cor- 
respondence with the 850 /im continuum emission of Davis 



et al. (1999) extending southward from SMMll. Another 



dust lane is extending roughly NE from SMM3 and co- 
incides quite well with a "bulge" of higher velocity gas 
visible in both the ^"^CO and C^^O centroid velocity maps 
(Fig. H). Here, as suggested by Davis et al. ( 1999[ ), the 
high- velocity gas and dust continuum emission are maybe 
showing the action of a wind powered by the SMM 2/3/4 
group. 



1*1 CnW - 



1*05110" - 




I1-I195Q.01 



00' 



Fig. 2. Overlay of SCUBA 850 /im gray-scale im- 
age (kindly provided by C.J. Davis, Davis et al. 1999) 
and C"'^*O(l-0) integrated emission (contour levels as in 
Fig. |l|). The thicker grey contour corresponds to 50% of 
the maximum. 



3.1.2. CS(2-1) and C34S(2-1) 

To determine the distribution of the denser gas in Serpens 
we obtained integrated intensity maps of the CS and 
C'^^S(2-1) emission (see also Sect. |3.4D . The maps are pre- 
sented in Fig. ^, where CS and C'^'*S(2-1) show several 
condensations and the two main clumps, NW and SE. A 
detailed comparison with the N2H+ and C^^O unveils sev- 
eral important differences. TSOO showed that the N2H+ 
map follows the spatial distribution of the submillimeter 
sources in both the NW and SE sub-clumps much more 
closely than the other molecular tracers such as C^^S and 
CH3OH. This is confirmed by our new single dish C'^^S 
and C^'^S maps, showing that the CS distribution resem- 
bles much more closely that of other shock tracers, such 



as H2CO (McMuUin et al. 200C), in agreement with what 
found by TSOO. 

3.2. Line profiles 

TSOO mapped the NW and SE clumps of the Serpens re- 
gion in the N2H"'"(l-0) transition and found four main 
cores of emission, A, B spatially associated with the 
SMM3/4 sources in the SE clump, and C, D associated 
with SMM1/S68N in the NW region. As mentioned ear- 
lier, they also identified an extended, lower surface bright- 
ness Spur, stretching from the NW to SE clumps. Spectra 
of N2H+ and also of CS, C^^S and C^^O at the emis- 
sion peaks of these regions are presented in Figs. ^ to |[ 
Our maps of the Serpens region obtained using CS, C'^^S 
and C^^O complement the N2H+ observations of TSOO 
and show some new features. In particular, the Spur is 
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Fig. 3. Maps of the CS and C^'^S integrated intensity 
in the velocity interval 2.0 to 13.0 km (CS) and 5.0 
to 11.0 km s-i (C^^S). In the CS map the value of the 
first contour is 0.24 K km s~^and the contour interval 
is 0.6 K km s^^. In the C^''S map the value of the first 
contour and the contour interval are both 0.09 K km s~^. 



detected in all molecular tracers observed by us and its 
emission only covers a narrow (~ 0.4 — 0.5 km s""'^) ve- 
locity range. The N2H+(l-0) emission is mostly optically 
thin, except at the emission peaks in the NW clump where 
r(N2H+) ^ 1. 



3.2.1. i'^CO(l-O) and C^^O(l-O) 

The ^■^CO spectra in the NW clump appear to be self- 
absorbed (see Fig. as the dip is at about the same 
velocity as that of the CS spectra (~ 8.6 km s~^) and is 
aligned with the C^'^S peak. The C^^O spectra appear to 
have a Gaussian profile in the NW clump and are mostly 
optically thin. In this clump the peak of the C^^O line is 
generally well aligned with those of C^'^S and N2H"'" and 
traces the dip of the ^^CO and CS spectra (see Figs. || and 
I)- 
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Fig. 4. Line spectra at position A(133,-177) and B(88,- 
44) in the SE sub-cluster. From top to bottom: CS(2-1), 
C34S(2-1), Ci«O(l-0), "CO(I-O), and N2H+(loi-0i2). 
The vertical lines are drawn at the peak velocities of the 
N2H+ Gaussian fits. The highest velocity N2II+ compo- 
nent corresponds to the emission from the "Spur" . 



In the SE clump the line shape of ^'^CO tends to be 
non-Gaussian, especially at offsets 6 < —100" and a > 20" 
from the map center. As for the N2H+ line shape, the 
C^^O non-Gaussian profile is due to the presence in the 
beam of several velocity components (cores "A" , "B" and 
the "Spur" of TSOO), as shown in Fig. ^ where the two 
N2II+ peaks are nicely aligned with those of C^^O. It is in- 
teresting to note that the ratio of ^^CO to C^^O emission 
is lower in the Spur than at the other positions. For exam- 
ple at the position and velocity of core A, (133,-177) and 
Vlsr = 7.15 km a'^, the "CO to C^O intensity ratio 
is about 3, whereas at the velocity of the Spur this ratio 
decreases to about 2. 

3.2.2. CS(2-1) and C34S(2-1) 

Most of CS spectra in the NW clump (see Fig. ||) are self- 
absorbed and are optically thick, as confirmed by the pres- 
ence of the C^'^S emission peaks at the dip of the CS spec- 



tra (see also Williams and Myers 1999). The differences in 
the profiles between the two isotopomers are mainly due to 
their different optical depths: the C'^'^S lines are optically 
thin and single-peaked whereas the CS spectra generally 
present two peaks. The asymmetries in the self-absorbed 
CS sp ectra have been interpreted by Williams and Myers 
( |1999|) as the signpost of a large-scale contraction of the 
gas around the clusters. The blue asymmetric line profile 
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Fig. 5. Same as Fig. || for positions C(-44,44) (i.e., 
SMMl) and D(-44,133) (or S68N) in the NW sub-cluster 
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Fig. 6. Same as Fig. | for positions E(-133,133) (or 
S68NW) and at one position in the Spur, (133,-310) 

of CS is visible in both the NW and SE clumps, and by 
comparing CS and C'^^S we will show later that in these 
regions the centroid velocity of the optically thicker iso- 
topomer is lower (bluer) than its rarer and optically thin 
counterpart. 



The CS spectra show line wings due to outflows from 
several cluster sources, such as S68N and SMMl in the 
NW clump, and SMM3 and SMM4 in the SE clump. The 
CS spectra in the Spur are self-absorbed, as shown in 
Fig. 1^, and both the N2H+ and C^^O emission peaks are 
aligned with the dip of the CS spectra. 

3.3. Kinematics of the Serpens region 

The analysis of the gas kinematics in the Serpens region 
is complicated by the presence of several velocity fields 
that may include infall, rotation, turbulence and outflow 
motions. In particular, the presence of several known out- 
flows may render the unambiguous detection of either ro- 
tation or infall very difficult. In this work we did not at- 
tempt to detect infall motions around some of the specific 
submillimeter sources, but we were instead interested in 
the global motions of the cloud and within the cloud. The 
molecular tracers used also reflect this goal: their exci- 
tation conditions make ^^CO and C^^O(l-O) less prone 
to contamination by outflows. Likewise, the two optically 
thin tracers C34S(2-1) and N2H+(l-0) are not affected by 
the presence of outflows. However, the CS(2-1) emission 
is known to be contaminated in the presence of strong 
outflows (see TSOO and refrences therein) and indeed we 
detect line wings in the CS lines (see Sect. 



3.2.2) 



We present the velocity information in four different 
graphical forms designed to emphasize different kincmat- 
ical properties of the gas: (i) line core, red and blue wing 
maps to show large-scale velocity structures; (ii) centroid 
velocity maps, to show small velocity trends of the cloud 
bulk emission; (iii) normalized centroid velocity difference 
between an optically thick and a thin tracer, to show sys- 
tematic differences such as infall motions; (iv) line width 
maps, to show the velocity dispersion accross the region 
for the various tracers. 



3.3.1. Line core and line wings maps 

To determine the velocity distribution of the ambient 
gas in Serpens we made integrated intensity maps of the 
blucshifted (5.0-7.3 km s'^), line core (7.3-8.4 km s^^) 
and rcdshifted (8.4 — 12.4 km s^^) emission of ^^CO and 
C^*O(l-0). The maps are shown in Fig. ^ where the 
^^CO and C^^O emission is clearly separated and occu- 
pies mainly the half-plane to the East in the blucshifted 
channel, whereas the emission is mostly concentrated in 
the half-plane to the West in the redshifted channel (see 
also Fig. H). In particular, ^^CO shows a very different 
distribution from that of CO (l-O) (McMuUin et al. |1994D 



and (2-1) (Davis et al. 1999). In fact, while the CO emis- 



sion is affected by the high- velocity gas radiated from the 
submillimeter sources, ^^CO is less affected by molecular 
onflows and thus can better show the kinematics of the 
large-scale gas. In the central velocity map of C^^O the 
Spur is clearly visible as a gas filament extending south- 
ward from the SE cluster of submillimeter sources. 
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The maps in Fig. ^ show emission from the SE and 
NW sub-clusters at all velocities, although the C'^^S emis- 
sion is weak in the 5.0 — 7.3 km s^^velocity interval and is 
mostly concentrated in the SE clump. We also note that 
the C^^O(l-O) emission from the NW clump is concen- 
trated in the line core and especially in the redshifted 
channel. The morphology shown by the channel maps of 
the two CO isotopomers and of the density sensitive trac- 
ers is clearly different. The ^^CO and C^*O(l~0) line emis- 
sion show a velocity gradient roughly E to W whereas 
CS(2-1) is tracing the embedded denser cloud cores and 
the EW velocity gradient is much less evident. 

Moreover, they show other isolated cores of emission, 
such as S68N, visible in the CS central velocity map. This 
core has a weak N2H+ counterpart shown in the N2H"'" 
channel map of TSOO at a velocity of about 8.2 km s~^. 
Another emission peak at the position of SMMl is visible 
in both the CS and C^'^S 8.4-12.4 km s^^channel maps; in 
C^^'S this emission peak is also visible in the other velocity 
channels. The Spur is also visible in the central velocity 
map. Finally, in the CS 8.4-12.4 km s^^channel map the 
SE sub-cluster appears to be composed of two separate 
cores, one at the position of SMM4 and the other slightly 
north of the SMM 3/6 group. 

The central velocity ^^CO map shows an emission peak 
at the position of SMM4 and of the CS peak. Another 
peak is visible sligthly NW of SMM3; the N2H+ core B 
is positioned between these two emission peaks. In the 
5.0-7.3 km s~^channel map a sharp peak of ^"^CO emis- 
sion is visible South of SMM4, which is also the region 
of C^^O emission, although the peaks are not coincident. 
A secondary peak of C^^O emission can be found in the 
NE region of this velocity channel that is coincident with 
enhanced ^'^CO emission. There are clearly two velocity 
components overlapping at this position, as shown by the 
C^^O and i^CO spectra. 

3.3.2. Large-scale velocity gradient 

The velocity gradient of an interstellar gas cloud can be 
determined by using all or most of the data in a map at 
once, by least-squares fitting maps of line-center velocity 
for the direction and magnitude of the best-fit velocity 
gradient (see Goodman et al. 1993 and references therein) . 

A cloud undergoing solid-body rotation would exhibit 
a linear gradient, dV/dr, across the face of a map, per- 
pendicular to the rotation axis. Thus, following Goodman 
et al. ( 1993 ) we fit the function 



Table 2. Results of velocity gradient fitting. 



Vlsr ^Vo + aAa + bA5 



(1) 



to the peak velocity of a Gaussian fit to the emission pro- 
file of several molecular tracers. Here, Aa and AS repre- 
sent offsets in right ascension and declination, expressed 
in radians, a and b are the projections of the gradient per 
radian on the a and 6 axes, respectively, and Vo is the 
LSR systemic velocity of the cloud. The magnitude of the 



Lino 


Vo 

[km s-^] 


dV/dr 
[km pc~^] 


0v 

[deg E of N] 


C^«O(l-0) 


8.31 


1.25 


-88.9 


^=^00(1-0) 


8.03 


0.75 


-70.8 


C^^S(2-1) 


8.42 


2.03 


-78.9 


N2H+(l23-0l2) 


8.34 


0.54 


-72.4 



velocity gradient, in a cloud at distance d, is thus given 

by 



dV _ ja^+b^y/^ 
dr d 

and its direction (measured east of north) is 

9v = tan~^ — 


In order to carry out a least-squares fit to Eq. (|l|) we have 
computed the following error function 



(2) 



(3) 



E 

hi 



VLSRihj) - Vo- aAui - bA6j 



(4) 



where Vi^snihj) is the velocity of the Gaussian fit to the 
spectrum at position (Aai, Adj) in the map, and ay{i,j) 
is the uncertainty in Vlsr determined by the fit, given by 
(see Goodman et al. 1993| and references therein): 



a„ = 1.15 ( ^) (<5„Ay)i/2 



(5) 



where T is the peak of a Gaussian fit to the line profile, 
ctt is the RMS noise in the spectrum, 6^ is the velocity 
resolution of the spectrum and is the FWHM line 
width. Only those points for which SNR= T/ar > 3 were 
used in Eq.^). 

To find the minimum of the function expressed by 
Eq.(^ we have used the dow nhill simplex method, using 
"AMOEBA" of Press et al. ( |1992D as the basic multidi- 
mensional minimization routine. The results are listed in 
Table || for the various molecular tracers. Two comments 
are in order: (i) N2H+ and C^'^S trace the denser gas in the 
cores and are thus less suited to calculate the large-scale 
velocity gradient; (ii) ^^CO is mostly optically thick with 
a non-Gaussian line profile and it is also not reliable. For 
these reasons we choose C^^O as the best available tracer 
to us of the velocity gradient in the Serpens cloud. Despite 
the aforementioned caveats, the direction of the velocity 
gradient is not significantly different from one tracer to 
another and confirms our early assumption in Sect. 3.3.1 
that the Serpens cloud exhibits a velocity gradient roughly 
E to W. 

3.3.3. Centroid velocity maps 

The centroid velocity of a line is that velocity at which 
the integrated intensity is equal on either side of the line 
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Fig. 7. From top to bottom: line core and line wings integrated intensity maps for C O, CO, C S, and C S. 



profile, and we compute it as in Narayanan et al. ( 19981 ). In Fig. ^ we present the centroid velocity maps of the 

Centroid velocity maps are usually a good indicator of ^^CO and C^^O lines. To minimize the potential contam- 

global velocity structures, such as rotational motions, and ination of the kinematics of the cloud by the outflow ve- 

we use them here as further evidence of the large-scale locity fields, the centroid velocities were computed over 

velocity gradient found in the previous section. line core intervals (6.0 — 10.5 km s~^) which must also be 
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wide enough to encompass the line core emission in both 
the NW and SE sub-clumps. Both maps show an increas- 
ing gradient in centroid velocity from roughly E to W, as 
expected. 

The maps shown in Fig. ^ also show a velocity "bulge" 
extending from about the map center towards the E, which 
may be caused by one of the redshifted CO(l-O) outflow 
lobes found by Narayanan et al. (2002), despite our cen- 
troid velocities are calculated using only line core emission. 
Even if some contamination of the centroid velocity maps 
by the outflows may still be present, the velocity gradient 
observed in the maps of Fig. || is extended to the whole 
Serpens cloud and is unlikely that the "'^^CO and C^^O 
centroid velocities are seriously affected by the outflows. 
We can also exclude that the observed velocity gradient is 
due to separate clumps with different systemic velocities, 
as the differences in the lines Vlsr are less or comparable 
with the line FWHM of C^^O and show a smooth change 
from E to W. As will be further discussed in Sect. 4.1, we 
conclude that the observed velocity gradient may indeed 
be caused by a global rotation of the Serpens molecular 
cloud whose rotation axis is roughly aligned in the SN 
direction. 

In Fig. H we also present the centroid velocity maps 
of the CS and C34S(2-1) transitions. The CS map shows 
essentially a plateau in the region where most of the CS(2- 
1) emission is detected. Our spatial resolution is not good 
enough to distinguish any particular feature around spe- 
cific submillimeter sources. However, we note once again 
the small bulge of redshifted emission E of the map center. 

3.3.4. Normalized velocity difference maps and infall 
speeds 

A diagnostics which is often used to disentangle the large- 
scale kinematics of molecular cloud cores is based on the 
normalized centroid velocity difference (see Mardones et 
al. 1997). The normalized centroid velocity difference SVc 



is a non-dimensional parameter defined as: 

K(CS)-K(C34S) 
'^'^ = AV(C3^S) 

and is very useful to quantify the observed line asymme- 
tries in the optically thick CS(2--1) line. Alternatively, one 
can use N2H+(l-0) to replace C^^S as the optically thin 
tracer in Eq. (||). A blue asymmetry {SVc < 0) can be 
interpreted as an indication of infall motion (Myers et 
al. |1996| ; Mardones et al. |1997| ). 

In Fig. ^ the map of SVc as defined in Eq. ^ is 
shown. We decided to use CS as the optically thick tracer 
and C^'^S, rather than N2H+, as the optically thin tracer 
to ensure that both species have homogeneous chemical 
properties. We can see that there are two regions of the 
map where the negative values of SVc are more concen- 
trated, which are approximately coincident with the SE 
and NW sub-clumps. This indicates that there is molecu- 
lar gas around the SE and NW sub-clusters which is still 
undergoing an infalling motion. 




Ac((arcsec) 

Fig. 8. Centroid velocity maps for: ^^CO, C^^O, C^^S, 
and C'^'^S (from top to bottom). The centroid velocities 
are obtained over the velocity interval corresponding to 
the line core (6 to 10.5 km s~^). Only the observed points 
where the integrated intensity is greater than three times 
its corresponding RMS uncertainty are used in calculat- 
ing the centroid velocity. The offsets are measured with 
respect to the center position listed in Table 0. 
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Fig. 9. Normalized centroid velocity difference map of 
CS and C'^'^S. The centroid velocities were calculated over 
the velocity interval 5.5 to 10.5 km s~^. To better show 
the contrast between the regions where SVc < and those 
where 5Vc > the upper limit for this map has been set 
equal to 0. The thick white line represents the contour for 
the 50% of the peak of the H2 column density. 



Typical infall speeds at different positions in the SE 
and NW clumps can be calculated by using the simple 
two layer model of Myers et al. (1996), which gives: 



red 



Vblu 



■log 



1 + cTbd/Tu 
1 + cTrd/Td 



(7) 



where is the brightness temperature of the dip in the 
CS line shape, Tbd is the height of the blue peak above 
the dip, and Trd is the height of the red peak above 
the dip. The velocity dispersion can be obtained from 
the FWHM of an optically thin line, such as N2H+, thus 
fJv — A^N^H+ZV^ log 2. The results are shown in Table ||, 
where the first two positions refer to the SE clump, and 
the next three correspond to the NW clump. In the SE 
clump (including position B, not shown in Table ^ the 
CS line profiles are more symmetric, implying a lower in- 
fall speed, except near position A where Vin > 0.2 km s~^. 
This suggests different states in the core contraction pro- 
cess, as also shown by Fig. |[ 

A similar difference clearly exists in the NW sub- 
clump, where positions C and E show substantially differ- 
ent infall speed compared to position D. This difference 
has also been observed by Williams & Myers (1999), and 
our results confirm their conclusion about the presence of 
supersonic (Mn > crth[H2] = 0.27 km s'^at T = 20 K) 
inward motions in the NW sub-clump, and indicate near- 
supersonic contraction in the SE sub-clump. 

For these same positions we can quantify the turbulent 
motion by comparing the N2H+ line FWHM with that 



Table 3. Line velocity dispersion of N2H+ (second 
column), non-thermal velocity dispersion (third 
column) and infall speeds (fourth column) towards 
the SE and NW sub-clumps (a = Ay/V81og2). 



Position 




(TNT 


Vin 




[km s~^] 


[km s^^] 


[km s"^] 


A 


0.33 


0.32 


0.26 


(88,-177) 


0.35 


0.34 


0.21 


C 


0.49 


0.48 


0.35 


D 


0.35 


0.34 


0.09 


E 


0.62 


0.58 


0.44 



expected for gas having equal thermal a nd no n-thermal 
motions (Mardones et al. 1997, Lee et al. 1999): 



AVo = 



;iog2fcr 



1/2 



(8) 



where T is the gas kinetic temperature, k is the Boltzmann 
constant, m^^ii+ is the molecular weight of N2H+ (29 
amu), and TOav is the gas mean molecular weight (2.3 
amu). For gas with T = 20 K, we find that AI442H+ > 
AVo = 0.66 km s^^, thus having greater non-thermal 
than thermal motions (see Table ^ . This result is consis- 
tent with the prevalence of turbulent infall motions in a 
sample of co res with embedded Y SOs fo und by Mardones 
et al. (1997), whereas Lee et al. ( |1999 ) found that most 
starless infall candidates are in a state of thermal infall. 
Incidentally, we note that also the starless core at position 
E has a greater turbulent rather than thermal motion. 

The aim of this paper is not to calculate the mass 
infall rate for every cloud core, as our spatial resolu- 
tion is not high enough. However, for consistency we did 
calculate the kinematic mass infall rate, (dM/di)kin = 
47ri??j {ma.^n-a^) Vin, toward positi on E. The infall radius 
was taken from Williams & Myers ( |1999| ), i?in = 6200 AU, 
Vin = 0.44 km s~^ from Table |[ and (mavfiHa) was ob- 
tained by taking the average of all mass estimates in 
Table |^ (except for the values obtained from the sta- 
tistical equilibrium method, see Sect. |X|) and divid- 
ing by the total volume of the NW clump. Thus, we 
obtain (dM/di)kin ^ 510"'^ Mgyr^^ This value can 
be compared with the gravitational mass accretion rate 
of (dM/dt)gr = 0.975a VG (Shu |l97^, where a is the 
velocity dispersion of the molecule of mean mass Wav, 
— kT/rriiiv + cr^-p. Using T = 20 K and ctnt from 



Table | we find (dM/di)gr ~ 5.9 10"^ Mgyr"!. The 
agreement of these two rates indicates that the derived 
inward motion at position E is consistent with gravita- 
tional infall. 



3.3.5. Line width maps 

Another useful diagnostic tool is to plot a measure of 
the width of the line profile. Given the deviation from 
a Gaussian line profile presented by some of the observed 



L. Olmi and L. Testi: The Serpens molecular cloud and protocluster 



11 



molecular tracers, in Fig. [l^ we show the maps of the 2"'* 
moment of the data. 

The map of the ^'^CO line width in Fig. |l^ clearly 
shows that the FWHM in the SE clump is enhanced com- 
pared to the rest of the ambient cloud. Although to a 
lesser extent, this is also seen in the NW clump. The C^^O 
line widths are also larger in the SE clump, which is also 
an effect of the presence of the two velocity components 
(clumps A, B and the Spur) in this part of the cloud. A 
two-component Gaussian fit in the SE region shows that 
most of the contribution to the line width of C^*O(l-0) 
comes from the lowest velocity component, although the 
peak of the line-FWHM is shifted relatively to the nominal 
positions of the clumps A and B. 

Likewise, the N2H+(l-0) map shows an enhacement of 
the line width in the SE clump, whereas this is less evident 
in the NW clump, also due to the limited extension of the 
map in that region. The CS(2-1) and C^''S(2-1) maps 
of the 2"'' moment of the data shown in Fig. |l^ show 
some enhacement of the line width at the positions of the 
submillimeter sources, but the most noticeable feature is 
the much larger line width at the position of the S68NW 
starless core described by Williams & Myers (1999). 



500 



3.4. Physical parameters 

In this section we derive the physical parameters of the 
molecular gas, using various tracers and approximations. 
We first derive estimates of the gas excitation and column 
densities using a simplified LTE approximation, then we 
derive particle densities from the C'^^S line intensities us- 
ing an LVG approximation combining our data with the 
C34S(5-4) observations of McMuUin et al. ( |1994| ). Finally, 
we provide various estimates of the masses of the entire 
cloud and of the various clumps identified by TSOO. 



3.4.1. Excitation and column densities 

The C^^O emission can be used to trace the column den- 
sity of the Serpens cloud. To derive values for the col- 
umn density we first estimated the excitation temperature 
throughout the entire region, based on the ^'^CO line tem- 
peratures at a velocity of 8.58 km s~^. We also estimated 
the C^^O optical depths using the ^^CO to C^^O line ra- 
tio and a relative abundance X[13qq]/j^[c18q]^ j^2. We 
then used both Tex and t(C^*0) to calculate the column 
density of a J ^ J — 1 rotational transition according to 
the following formula (Lis & Goldsmith 1991 ): 



4.0 X IQi^ fEj 



1 



r?mb 1 



Tldv (9) 



where B denotes the rotational constant, Ej is the upper 
state energy, fj, is the dipole moment and we use the es- 
cape probability t/[1 — exp(— r)] to account for first order 
optical depth effects. We also assume a unity beam-filling 
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Fig. 10. Maps of the 2'"^ moment of the data (line width). 
From top to bottom: ^^CO, C^^O, N2H+, C^^S, C^^S. 
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factor. The excitation temperature is obtained by invert- 
ing the equation of transport: 



Tnib = [J. (Tex) - J.(Tbg)] [1 - exp(-r)] 



(10) 



where Tex and r are calculated at a specific velocity and 
T is evaluated as described above. Moreover, for fcTex >> 
hB the partition function of a linear molecule is given by 



hB 



(11) 



where h is the Planck constant. The distribution of 
-^tot(H2) calculated using C^*0 data is shown in Fig. pT| 
and corresponds very well to the map of the C^^O inte- 
grated emission of Fig. |l|. 

The CS and C'^^S data can be used to derive the gas to- 
tal column density following the method described above 
for ^^CO and C^^O. C'^'^S is in fact optically thin and can 
thus be successfully used to trace the column density of 
the gas if its relative abundance can be determined. We 
derived the C'^^S abundance using the LVG approxima- 
tion to model the measured line intensities, as described 
later in this section. Therefore, we can use C^^O and C'^'^S 
to compare the spatial distributions of the column den- 
sity and excitation temperature at two velocities char- 
acterizing the NW (8.58 km s'^) and SE (7.3 km s^^) 
sub-clumps. We find that CS and C^^O show a different 
behaviour: in fact, the C^^O maps present a spatial offset 
between the peaks of Tex and those of the gas column den- 
sity distributions (see Fig. ^l]). At the velocity of the NW 
sub-clump the C^^O column density map shows a clump 
that surrounds the NW group of submillimeter sources, 
whereas the map of Tex shows a smooth plateau at the 
same position. At the velocity characterizing the SE sub- 
clump these differences are less extreme: a broad peak of 
column density surrounds sources SMM2/4/11, whereas 
the Tex map peaks between SMM4 and SMM3/6. In both 
the NW and SE sub-clumps the distribution of the col- 
umn density (see Fig. ^ follows quite closely the map of 
the C^^O integrated emission shown in the bottom panel 
of Fig. 0, suggesting that the emission peaks are column 
density maxima and are not due to peaks in the excitation 
of the molecular gas. 

The CS and C'^'^S transitions can also be used to deter- 
mine the spatial distribution of Tex(C'^^S) and the column 
density. The map of Tex(C^^S) has a broad peak including 
the SMM 3/4/6/8 group. The gas traced by CS, however, 
is characterized by spatial distributions of Tex and column 
density that are much more similar (see Fig. |l^). In the 
NW sub-cluster the column density map shows a peak at 
the position of S68NW, while Tex has a maximum slightly 
NW of the S68NW position. At the velocity characteriz- 
ing the SE sub-cluster, the Tex and column density maps 
both have a peak very close to the position of SMM4 and 
show little differences otherwise. Therefore, there is a ten- 
dency of the gas column density to have a different spatial 
distribution from that of Tex in the NW sub-cluster, and 
this difference is much more evident in C^^O. 



C O V=7.3 km/s 




-500 500 
C'^O V-8.58 km/s 




500 500 



Aa(arcsec) 



Fig. 11. Maps of Tex and of H2 column density 
using C^^O and a velocity of 7.3 km s^^(top) and 
8.58 km s~^ (bottom), characteristic of the SE and NW 
clumps, respectively. 



3.4.2. Density estimates 

The physical number densities towards the densest con- 
densations can be derived using density tracers such as the 
optically thin C^'*S transitions. Because we observed the 
J = 2 - 1 line only, we also used the C3'^S(2-1) and ( 5-4) 
observations towards SMMl of McMulhn et al. ( |1994D . By 
using the J = 2 — 1 measurements taken with the FCRAO 
and NRAO telescopes we were able to estimate the decon- 
volved angular diameter of the source 9s — 70". Assuming 
a gas temperature of 35 K towards SMMl and a velocity 



gradient of about 21 km s~ pc~ , we then used the LVG 
approximation to model the measured line intensities, ob- 
taining a number density rifja — 5.6 x 10^ cm~^ and a rela- 
tive abundance a: [C^^S]/ a: [H2]= 4x 10"", which are con- 



sistent with the values obtained by McMullin et al. (2000). 
There are no C^'*S(5-4) observations towards the SE sub- 
clump. However, an estimate of the physical conditions 
at the SMM4 position can be derived by using the LVG 
model parameters obtained for SMMl, and changing one 
parameter only, e.g. the density, to fit the C^''S(2-1) line 
intensity. By doing this we derive — 7 x 10^ cm^'^ 
towards SMM4. 



3.4.3. Mass estimates 

Integrating over an area of 111 arcmin^, or 0.90 pc^, the 
total mass of the C-'^^O emitting region is about 200- 
500 Mq depending on the adopted value for the excitation 
temperature and assuming a distance to Serpens of 310 pc 
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where Av is the average full width at half maximum of 
the line (-1.3 and ~1.6 km s-^for C^^O and C^^S, re- 
spectively). The typical diameter of the source is obtained 
as = where fig is the solid angle, within the 50% 

contour, subtended by each sub-clump and measured us- 
ing Figs, [l] and |l|. The masses of the two condensations 
can also be constrained assuming a spherical geometry 
and the peak density derived from a statistical equilib- 
rium model applied to the C'^'^S transitions (see above), 
and using 

MsaiMe) = 2. 926x10"^ d^(kpc)6',^(arcsec)nH2(cm-3)(14) 

Clearly, Mgc should be considered an upper limit to the 
mass of the molecular gas as we have used the peak and 
not average density values, and it is also much more sen- 
sitive to the error on the distance to the cloud. In Table |^ 
we present a summary of the various mass determinations 
for the two sub-clumps SE and NW, including the values 
determined in TSOO and TS98 for the N2H+(l-0) virial 
mass and for the total mass of the dust cores and stellar 
components within the clumps. 



Fig. 12. Maps of Tex and of H2 column density 4_ Discussion 



using CS and a velocity of 7.3 km s^ (top) and 
8.58 km s^^ (bottom), characteristic of the SE and NW 
clumps, respectively. 



and a C^^O fractional abundance of X[Ci^0]/X[H2] = 
10-"^ (McMuUin et al. ^000| ). The ^^CO excitation tem- 
perature used to calculate the total mass i n the Serpens 
cloud can be estimated as described in Sect. 3.4.1 , obtain- 
ing values in the range 



4.1. Large-scale cloud kinematics 



As discussed in Sect. 3.3, our ^■^CO and C^^O data show 



a large-scale gradient in the line centroid velocity, which 
appears to be consistent with what expected if the cloud 
is in slow rotation. Using the C^*0 value for the velocity 
gradient found in Sect. 3.3.2 we find an angular velocity 
w ~ 4.1 X 10-^'^ s-i and following Goodman et al. (|1993D 



_ 7 to 22 K, used to calculate the 
column densities in Fig. ^ However, these temperatures 
are lower than the gas temperatures derived by McMuUin 
et al. (2000), due to the fact that the peak ^^CO line tem- 
perature traces the cooler larger scale cloud. The range in 
the cloud total mass given above reflects these uncertain- 
ties in the value of the excitation temperature. 

The maps shown in Figs. ^ and [l^ are likely to un- 
derestimate the column density when the ^^CO excitation 
temperature is used. Therefore, we have estimated the col- 
umn density using an overall gas temperature of 25 K from 
McMuUin et al. (200C), it is then possible to obtain the 
masses of the molecular gas in the NW and SE clumps 
using the formula: 



we can then calculate the parameter /3, defined as the ratio 
of rotational kinetic energy to gravitational energy. For a 
sphere with constant density, po, (3 can be written as: 

wehere G is the gravitational constant. To calculate po 
we use the total gas mass, ^--^300 M0 (corresponding to 
an extended region that includes both the NW and SE 
clumps and the lower density envelope) and the corre- 
sponding average radius, 0.55 pc. Hence, po — 2.9 x 
10~2° g cm-3. Substituting in Eq. ( pi) we find /3 ~ 0.07. 



Afed = m(H2) Mg / A^(H2)dr! 



(12) 



where / N{)l2) dll is the hydrogen column density inte- 
grated over the region enclosed by the contour level at 
50% of the peak value for each sub-clump, /ig is the mean 
molecular weight (1.36), r7i(H2) is the mass of a hydrogen 
molecule and d is the distance to the Serpens cloud. The 
masses derived in this way (see Table ^) can be compared 
with the virial masses, 

M„(Mq) = 0.509 d(kpc) 6',(arcsec) Av^{\tcl &-^) (13) 



All clouds studied by Goodman et al. (1993) and Barranco 
& Goodman ( 1998| ) have j3 < 0.18, which is low enough 
to inhibit, according to those authors, fragmentation due 
to rotationally driven instabilities. Hence, the fragmenta- 
tion observed in Serpens may have involved both gravi- 
tational and magnetic instabilities, but it is very unlikely 
that the cloud rotation could have induced fragmentation. 
Incidentally, we note that the direction of the velocity 
gradient is approximately perpendicular to the average 
outflow axis determined by TSOO (see also Table |). It 
is tempting to speculate that the global cloud rotation 
may have influenced the orientation of circumstellar disks 
surrounding each individual (proto-)cluster member; how- 
ever, at this time this is a pure speculation since we can- 
not prove whether there is a physical connection between 
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Table 4. Med is the mass obtained from the column density integrated over an area corresponding to 50% of the 
peak value. M^ii is the virial mass, also obtained using the diameter of the 50% contour level, and Mge is the mass 
of a spherical clump with uniform density equal to the peak hydrogen density determined from statistical equilibrium 
models. Md is the total mass of compact dust cores and is the estimated total mass of near infrared stars. 



Sub-clump 






cd 
















(Mo) 


(Mo) 


(Mq) 


(Mq) 


(Mo) 


(arcsec) 


(Mo) 


(Mo) 




NW 


60 


38 


54 


58 


130 


139 


27 


19.5 


4 


SE 


50 


41 


80 


65 


220 


154 


~30 


~12 




Assume To, 


= n = 


= 25 K (McMuUin et al. 


200C|), and 


Aw ~1.C 


) kms"^ 






'''' Assume To 


= rk = 


= 30 K (McMuUin et al. 


2000D, and 


Av ~1.( 


) kms~^ 







t'^' Values from TS98, TSOO, and references therein 



Table 5. Outflows directions position angles from the 
compilation of TSOO 



Source 


Outflow direction 
[deg E of N] 


Outflow tracer 


Reference^""' 


S68N 


140 


CS, CH3OH 


1,2 


A3 


140 


Ha 


3 


SMM5 


150 


Refl. Nebula 


4 


SMMl 


140 


Ha, CS, CO 


1,2,3 


SMM3 


170 


Ha, CS 


1,5 


SMM4 


180 


Ha, CS 


1,6,7 


Average 


~155 







'•"^ 1: TS OO; 2: Wolf-Chase et a l (|l998|) ; 3: Hodapp ( |l999| ) 
4: Kaas (|l999| ); 5: Herbst et al ( |1997| ); 6: Hogerheijde et 
al. ([1999); 7: Eiroa et al. (|l997|). 



these two similar directions. The relationship between the 
large-scale velocity gradients and the young stellar objects 
embedded in the Serpens core should be further studied by 
measuring the velocity gradients in the circumstellar disks 
surrounding most of the objects (see e.g. Hogerheijde et 
al. 



1999) 



The simple rigid body rotation approximation used 
to derive to and f3 may be questioned in a cloud core 
that displays a large non-thermal contribution to the line 
width, an indication of complex motion on scales much 
smaller than the beam width. Nevertheless, Burkert & 
Bodenheimer ( 2000] ) have shown that turbulent cloud 
cores with random gaussian velocity fields may have a net 
angular momentum different from zero and may display 
a global velocity gradient. Surprisingly, the values of the 
specific angular momentum derived assuming a rigid body 
rotation represent a good approximation to the real values 
induced by the turbulent motions, at lea st within a factor 
of a few (Burkert & Bodenheimer 2000| ). Our conclusion 
that the measured cloud angular momentum is not large 
enough to play a dominant role in the fragmentation of the 
cloud and on the formation of the cluster is thus robust. 

At present we cannot exclude that the observed veloc- 
ity gradient is produced by a more complex kinematical 
structure of the cloud, such as a shear motion or a non- 
spherically symmetric contraction of the cloud. A shear 



motion would imply a net angular momentum of the cloud, 
similar to the rigid rotation discussed above. In fact, a 
shear motion can be seen as a complex rotational pattern, 
and would not change significantly the discussion above. 

Perhaps more interesting would be the possibility of 
a global, non sperically symmetric cloud contraction. A 
comparison of the total cloud mass as derived from the 
C^^O column density (^300 Mq) and the virial mass 
(~150 Mq), suggests that a global contraction may be 
possible. This would also be consistent with the infall sig- 
nature detected in the innermost regions of the cloud (see 
Sect. |3.3| and ^4[ ). However, the measured velocity gradi- 
ent (see Table ^ and Fig. ||) would imply peak contraction 
velocities <; 0.5 kms~^, which is larger than the infall ve- 
locities measured in the inner regions of the cloud (see 
Table ||) . We thus conclude that this is an unlikely possi- 
bility. 

4.2. Cloud contraction and large-scale physical 
conditions 



In Sect. 3.4 we showed that the cloud core displays two 
column density maxima roughly coincident with the two 
sub-clump identified by TSOO. There is also a strong 
correlation between the high column density regions and 
the infalling signature in the normalized centroid veloc- 
ity difference map (Sect. 3.3 and Fig. ||). As noted in the 
higher resolution N2H+(l-0)/CS( 2-l) study of the NW 
sub-clump by Williams & Myers ( 2000| ), it appears that 
the infalling speeds are higher in the high column den- 
sity regions surrounding the already formed young stellar 
objects but not at their precise location, where the feed- 
back produced by protostellar outflows is acting against 
the infall. Thus, our observations, together with those of 
Williams & Myers, support the view that the conditions 
for collapse are met in the higher density regions of the 
cloud, which are globally contracting. At a few selected lo- 
cations the collapse has already produced YSOs that, by 
means of their molecular outflows, are starting to revert 
the accretion and disperse the cloud core. 

We have also calculated the turbulent crossing time for 
the main sub-clumps in Serpens. The crossing time is de- 
fined as the radius divided by the Gaussian dispersion in 
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internal velocity (see Sect. 3.3.4). We used the solid angle 
subtended by the 50% contour (of the corresponding peak 
value) to calculate the radius of the main N2H+ cores of 
emission, B, C and D. We repeated the same operation 
for the core E and for the SE sub-clump mapped in the 
C^'*S(2-1) line. The crossing times are remarkably similar, 
varying in the interval 1.2 x 10^ to 1.5 x 10^ yr, although 
the radii vary by a factor of 2 (0.055 pc to 0.115 pc) and 
the velocity dispersions by almost a factor of three (0.35 
to 0.98 km s~^). In the case of the NW sub-clump these 
values are consistent with the time needed to reach chem- 
ical steady state, i.e . seve ral times 10^ years according 
to McMuUin et al. (^GGO|) . These large cores coexist in 
Serpens with smaller and younger regi ons ( e.g., Williams 
& Myers |l999| , |2000| , Narayanan et al. |200^ suggesting a 
wide range of evolutionary stages. 

4.3. Star formation efficiency in Serpens 

Combining our data with the results of TS98, TSOO, 
Giovannetti et al. ( |1998D , Kaas ( |1999[ ), and White et 



al. (1995) we can estimate the average and local star for- 
mation efhciency in the Serpens cloud core. In Table | we 
summarize the mass estimates for the gaseous sub-clumps, 
for the dust cores and for the young stellar objects embed- 
ded within them. The total mass of the Serpens cloud has 



been estimated to be in the ra nge 300 (see Sections ^T 



3.4.3 and McMuUin et al. gOOOj) to 1500 Mq (White et al 



1995 ), whereas the total mass of the young stellar cluster 
is estimated to be ^^25-40 Mq to which one should add the 
mass of forming stars in the dust cores of TS98. Assuming 
that the cores will produce stars with a 50% efhciency, 
meaning that half of each core mass will end up in a young 
star and half will be ejected in a protostellar outflow, the 
total mass of the final cluster is estimated in the range 
40-60 Mq. The global star formation efficiency is thus es- 
timated in the range 3 to 10%. However, within the two 
sub-clumps defined by the N2lI+(l-0) FWHM contours, 
where the gas column density is higher, the global infall 
motion is detected (Sect. ^ and Wilhams & Myers ^OOOj ) 
and where most of the submillimeter prestellar cores and 
protostars are located (TS98, TSOO), the Zoca/ star forma- 
tion efficiency is much higher (25 to 50%), suggesting that 
where the collapse conditions are met the star formation 
process proceeds with a high efficiency. 

5. Conclusions 

We performed wide-field imaging (up to about 16' x 16') 
of the emission towards the Serpens cloud core using mod- 
erately optically thick (13CO(1-0) and CS(2-1)) and op- 
tically thin tracers (Ci8O(l-0), C34S(2-1), and N2H+(1- 
0)). Our main goal was to study the large-scale distribu- 
tion of the molecular gas in the Serpens region and to 
understand its relation with the denser gas in the cloud 
cores, previously studied at high angular resolution. 

The ^^CO and C^*O(l-0) emission traces the distribu- 
tion of the column density in the Serpens cloud. The ^■^CO 



emission is extended and encompasses both the NW and 
SE clusters of submillimeter continuum sources, whereas 
the C^^O emission shows the clear demarcation between 
the NW and SE sub-clumps, which are also separated in 
velocity (see also TSOO). Maps of the OS and C^'^S{2- 
1) emission show several condensations in the two main 
sub-clumps, NW and SE, but present otherwise a differ- 
ent appearance from N2II+. The latter follows the spa- 
tial distribution of the submillimeter sources in both the 
NW and SE sub-clumps much more closely than the other 
molecular tracers. In both the NW and SE sub-clumps the 
distribution of the column density follows quite closely the 
map of the C^'^O and C'^'^S integrated emission, suggest- 
ing that the emission peaks are column density maxima 
and are not due to peaks in the excitation of the molecu- 
lar gas. We also find that the SE and NW sub-clumps are 
virialised with masses typically of order 60 Mq and peak 
number densities of order — 7 x 10^ cm^^^. 

The "CO and C"O(l-0) maps of the centroid veloc- 
ity show an increasing but smooth gradient in velocity 
from E to W. Using most of the data in a map at once, by 
least-squares fitting maps of line-center velocity for the di- 
rection and magnitude of the best-fit velocity gradient, we 
find that C^^O has a 1.25 km s^^ pc^^ velocity gradient 
which is directed from E to W. The velocity gradient is ex- 
tended to the whole Serpens cloud and is unlikely that the 
^^CO and C^^O centroid velocities are seriously affected 
by the outflows. We can also reasonably exclude that the 
observed velocity gradient is due to separate clumps with 
different systemic velocities, or that it is due to a global, 
non-symmetric cloud contraction. We thus think that the 
observed velocity gradient may indeed be caused by a 
global rotation of the Serpens molecular cloud whose ro- 
tation axis is very close to the SN direction, similar to the 
averag e outfl ow axis found by TSOO. Following Goodman 
et al. ( 1993 ) we also calculated the parameter (3, defined 
as the ratio of rotational kinetic energy to gravitational 
energy . We find (3 ~ 0.07 which, according to Goodman 
et al. ( 1993 ), is low enough to inhibit fragmentation due 
to rotationally driven instabilities. The cloud angular mo- 
mentum is not sufficient for being dynamically important 
in the global evolution of the cluster. Nevertheless, the 
fact that the outflows are roughly aligned with it may 
suggest a link between the large-scale angular momentum 
and the circumstellar disks around individual protostars 
in the cluster. However, our data does not allow to prove 
this connection. 

The normalized centroid velocity difference, SVc, is a 
typical indicator of underlying infall. Using CS and C^^S 
we find two large regions of the map where the negative 
values of SVc are more concentrated, which are approxi- 
mately coincident with the SE and NW subOclumps. We 
find that the SE and NW sub-clumps have greater non- 
thermal than thermal motions, and we calculate infall ve- 
locities of about 0.1 to 0.4 km s~^. Together with the 



results of Wilhams & Myers ( |2000|) , our observations sug- 
gest that infalling conditions are met in a relatively large 
region where the column density of the gas exceeds a cer- 
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tain threshold; when a protostellar object is formed, its 
outflow locaUy acts against the infall. 

Whithin the high column density regions, where the 
infalling signature is detected, the local star formation ef- 
ficiency is much larger than in the entire cfoud. This is 
consistent with the picture that where the conditions for 
star formation are met, then the process is highly effi- 
cient. Given that outflows from young protostars appear 
to oppose locally to the infall, such high star formation 
efficiency can be achieved only if most of the objects in 
the cluster are formed over a very short timescale. This 
short formation timescale, which does not allow much time 
for interaction between protostars given the relatively low 
density of objects observed in Serpens, and the indica- 
tion that the prestellar cores mass function is consistent 
with the field stars IMF (TS98), suggest that dynami- 
cal interactions and competitive accretion may not be the 
dominating processes in these types of regions. 

The observational connection of the large-scale rota- 
tion and the individual stellar systems flows, and the ex- 
tended slow inflow motion limited to the high column den- 
sity regions, appear to be consistent with a slow global evo- 
lution of the cloud that locally reachs the conditions for 
gravitational collapse and star formation. The very short 
timescale for the formation of the bulk of the cluster is 
advocated by bot h the accelerated star formation view of 
Palla & Stabler ( ^000 ) or the dynamical star formation 
view of Elmegreen ( |200C ) . 

It is difficult, from our observations, to give a definitive 
answer to whether star formation occurs in a slow quasi 
static fashion up to the point of instability or whether 
the entire process is driven by the fast dissipation of tur- 
bulence in clouds. Our data appear to be in qualitative 
agreement with the expectation of a slow contraction fol- 
lowed by a rapid and highly efficient star formation phase 
in localized high density regions. Our evidence is not 
conclusive, but we provide a series of observational con- 
straints with which numerical simulation have to be tested 
against quantitatively. From the observational point of 
view, on the one hand the possible connection between 
the large-scale rotation and the individual circumstellar 
regions should be made more precise by means of high 
angular observations, on the other hand comparisons be- 
tween wide area studies of the physical and kinematical 
conditions within giant molecular clouds and the distri- 
bution and properties of the young (forming) stars within 
them should be extended to other complexes. 
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